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Skyrme-QRPA calculations for low-lying excitation modes in 
deformed neutron-rich nuclei 
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Abstract. Low-frequency modes of excitation in deformed neutron-rich nuclei are studied by means of 
the quasiparticle random-phase approximation on the Skyrme-Hartree-Fock-Bogoliubov mean field. We 
investigate the microscopic structure of the soft = 0"*" modes systematically in neutron-rich Magnesium 
isotopes with A'' — 22, 24, 26 and 28 close to the drip line, and it is found that the strong collectivity in 
•^■^Mg and ''"Mg is acquired due to the coherent coupling between the f3 vibration and the pairing vibration 
of neutrons. Microscopic structure of the K'^ = 2"*" modes changes gradually associated with the location 
of the Fermi level of neutrons, and it is found that the proton particle-hole excitation generating the 
7— vibrational mode in '^''Mg continues to play a key role in the near-drip-line nucleus *°Mg. The low- 
frequency octupole excitations are also investigated and the microscopic mechanism for the enhancement 
of transition strengths is discussed. 



PACS. 21.60.Jz Nuclear Density Functional Theory and extensions - 
Collective models - 27.30. -l-t - 27.40. +z 



21. 10. Re Collective levels - 21.60.Ev 



1 Introduction 



Collective motion in unstable nuclei has raised a consider- 
■ able interest both experimentally and theoretically. This is 
I because low-frequency modes of excitation are quite sen- 
. sitive to the shell structure near the Fermi level, and we 
' can expect unique excitation modes to emerge associated 
with the new spatial structures such as neutron skins and 
the novel shell structures that generate new regions of de- 
formation Ij. 

In order to investigate new kinds of excitation modes 
in exotic nuclei, the random-phase approximation (RPA) 
based on the self-consistent mean field has been employed 
by many groups. (See Refs. ,2,3,4J for extensive lists of 
references concerning the self-consistent RPA and mean- 
field calculations. ) They are however largely restricted to 
spherical systems, and the low-frequency excitation modes 
in deformed neutron-rich nuclei remain mostly unexplored. 

Recently, low-lying RPA modes in deformed neutron- 
rich nuclei have been investigated by several groups [5ll6l[7l 
[51[51 [TUlllllll2lFra] . These calculations, however, do not take 
into account the pairing correlations, or rely on the BCS 
approximation for pairing (except for Ref. [13]), which 
is inappropriate for describing the pairing correlations in 
drip line nuclei due to the unphysical nucleon gas prob- 
lem |14| . Quite recently, we have developed a new frame- 
work of the self-consistent deformed quasiparticle-RPA 
(QRPA) based on the Skyrme-Hartree-Fock-Bogoliubov 
(HFB) mean field [H]. 



Presently, small excitation energies of the first 2+ state 
and striking enhancements of B{E2] 0^ 2^) in ^^Mg [TBI 
[T7] and ^''Mg [TSIITM^ are under lively discussions in 
connection with the onset of quadrupole deformation, the 
breaking of A'^ = 20 spherical magic number, the pairing 
correlation and the continuum coupling effects [^TE^ll^ 
[241125] . In order to get clear understanding of the nature 
of quadrupole deformation and pairing correlations, it is 
strongly desired to explore, both experimentally and the- 
oretically, excitation spectra of these nuclei toward a drip 
line [1117111111131311132] ■ 

In this paper, we apply the new calculation scheme to 
the low-frequency excitation modes in neutron-rich Mag- 
nesium isotopes close to the drip line, and investigate the 
microscopic mechanism of the excitation modes uniquely 
appearing in deformed neutron-rich nuclei. 

The paper is organized as follows: In Sec. [H the de- 
formed Skyrme-HFB -I- QRPA method is recapitulated. In 
Sec.jSl results of numerical analysis of the low- lying excita- 
tion modes in deformed neutron-rich Magnesium isotopes 
are presented. Finally, summary is given in Sec. |4l 



2 Method 

2.1 Skyrme-HFB in cylindrical coordinates 

In order to describe simultaneously the nuclear deforma- 
tion and the pairing correlations including the unbound 
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quasiparticle states, we solve the HFB equations [T4ll33l 

in coordinate space using cylindrical coordinates r = (p, z,(f 
We assume axial and reflection symmetries. Here, q = v 
(neutron) or tt (proton). For the mean-field Hamiltonian h, 
we employ the SkM* interaction [34, . Details for express- 
ing the densities and currents in the cylindrical coordinate 
representation can be found in Ref . |35| . The pairing field 
is treated by using the density-dependent contact interac- 
tion [551. 



5{r^r'). (2) 



where Qo{r) denotes the isoscalar density of the ground 
state and the spin exchange operator. Assuming time- 
reversal symmetry and reflection symmetry with respect 
to the x — y plane, we have to solve for positive Q and pos- 
itive z only, J7 being the z— component of the angular mo- 
mentum j. We use the lattice mesh size Ap = Az — 0.6 fm 



and a box boundary condition at p-^ 



9.9 fm, z-a 



12 



fm. The differential operators are represented by use of the 
11-point formula of Finite Difference Method. Because the 
parity and fl are good quantum numbers in the present 
calculation scheme, we have only to diagonalize the HFB 
Hamiltonian ([1]) for each 17'^ sector. The quasiparticle en- 
ergy is cut off at i^qp.cut = 60 MeV and the quasiparticle 
states up to fl'^ — 15/2^ are included. 

The pairing strength parameter tg is determined so 
as to reproduce the experimental pairing gap of "^"^Mg 
(Ziexp = 1.7 MeV) obtained by the three-point formula [37j . 
The strength t'^ = — 295MeV-fm'^ for the mixed-type in- 
teraction (tg = — IS.TSfg) [38] with j — 1 leads to the 
pairing gap {A„) = 1.71 MeV in ^'^Mg. 



2.2 Quasiparticle-basis QRPA 

Using the quasiparticle basis obtained as the self-consistent 
solution of the HFB equations ([T]), we solve the QRPA 
equation in the matrix formulation [39J 



E 

7(5 



fn 

9^.s 



1 

-1 



Jap 
9ap 



(3) 



The residual interaction in the particle-hole (p-h) chan- 
nel appearing in the QRPA matrices A and B is derived 
from the Skyrme density functional. We neglect the spin- 
orbit interaction term C^"^ as well as the Coulomb inter- 
action to reduce the computing time. We also drop the 
so-called " J^" term both in the HFB and QRPA cal- 
culations. The residual interaction in the particle-particle 
(p-p) channel is derived from the pairing functional con- 
structed with the density-dependent contact interaction 



Because the full self-consistency between the static 
mean-field calculation and the dynamical calculation is 
broken by the above neglected terms, we renormalize the 
residual interaction in the p-h channel by an overall factor 
/ph to get the spurious K'^ = 0~ and 1~ modes (repre- 
senting the center-of-mass motion), and I-C^ = 1+ mode 
(representing the rotational motion in deformed nuclei) at 
J zero energy (uph /ph'i'ph)- We cut the two-quasiparticle 
(2qp) space at Ea + Ep < 60 MeV due to the excessively 
demanding computer memory size and computing time for 
the model space consistent with that adopted in the HFB 
calculation; 2ii^qp^cut = 120 MeV. Accordingly, we need an- 
other factor /pp for the p-p channel. We determine this fac- 
tor such that the spurious K'^ — 0+ mode associated with 
the particle number fluctuation (representing the pairing 



rotational mode) appears at zero energy {vpp — *■ /, 



pp "^pp 



(See Ref. 15j for details of determination of the normal- 
ization factors.) 

In the present calculation, the dimension of the QRPA 
matrix ([3]) for the quadrupole K"^ = 0+ excitation in ^"Mg 
is about 16 400, and the memory size is 24.2 GB. The 
normalization factors are /ph = 1.106, and /pp — 1.219. 

In terms of the nucleon annihilation and creation oper- 
ators in the coordinate representation, ^{ra) and ^^p^ra), 
the quadrupole operator is represented as 

Q2K = J2 j drr''Y2K{f )^Hra)i,{ra). (4) 



The intrinsic matrix elements (?T.|Q2/f |0) of the quadrupole 
operator between the excited state \n) and the ground 
state |0) are given by 

{n\Q2K\0) = Q2kIpU2p + 9lp) = E <kIp- (5) 

a/3 aP 

The explicit expression of Q^2Kap given in Ref. [32] . The 
neutron (proton) matrix element M^, {M-^) is defined 

M. = M27lp, = <kIp- (6) 



ap£i' 



ap^TT 



2.3 Elimination of the spurious center-of-mass modes 

It is known that the self-consistent RPA, if the same effec- 
tive interaction or the same energy density functional is 
used exactly both for the ground state and for the excited 
state, restores translational invariance [40]. Because the 
present calculation scheme is not fully self-consistent, the 
calculated states \n) for — 0~ and 1~ excitations may 
contain the spurious component of the center-of-mass mo- 
tion. In order to separate the intrinsic excitations from the 
spurious excitation, the "physical" states \n) are assumed 



Xn = X^- XpP - XrR, 



(7) 



where Xn is an annihilation operator of the calculated 
RPA mode, R and P are the coordinate and momentum 
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operators of the whole nucleus {[R,P] — ih). The coeffi- 
cients X considered to be small because the spurious 
component is expected to reasonably decouple from the 
physical solutions. The physical states satisfy the follow- 
ing conditions. 

1. The vacuum condition: 

Xn\0) = (8) 

2. The decoupling condition: 

(0|[Xft,P]|0) =0,(0|[Xs,i?]|0) =0 (9) 
These conditions determine the coefficients Xp, X^', 

Xl = l{0\ [Xn, R]\0)/K Xr = -m[Xn, P] \0)/h. (10) 

The orthonormality of the physical states \h) satisfies in 
first order of the correction coefficient x; 



\[Xn,Xi]\0) 

\[x„,xl]\o)-zh{x-pxT-xWn 



(11) 



In the actual calculations, the correction coefficient 
is at most of order 10~^. Separation of the spurious modes 
is also proposed in Ref. [41] in a similar way to the tran- 
sition density. 

The explicit expressions for calculating the matrix el- 
ements of the octupole transition operator are given in 
Appendix. 



3 Results and discussion 
3.1 Ground state properties 

In Table [H the ground state properties are summarized. 
The neutron-rich Magnesium isotopes under investigation 



Table 1. Ground state properties of 34.36,38,40]yjg obtained by 
the deformed HFB calculation with the SkM* interaction and 
the mixed-type pairing interaction. Chemical potentials, de- 
formation parameters, average pairing gaps, root-mean-square 
radii for neutrons and protons are listed. The average pairing 
gaps of protons are zero in these isotopes. The average pairing 
gap is defined (zi), = — / drhg/ J drg. 





^-^Mg 


36 Mg 


38 Mg 


«Mg 


A„ (MeV) 


-4.16 


-3.24 


-2.41 


-1.56 


K (MeV) 


-19.8 


-21.0 


-23.7 


-24.4 




0.35 


0.31 


0.29 


0.28 




0.41 


0.39 


0.38 


0.36 


{A), (MeV) 


1.71 


1.71 


1.64 


1.49 


^{r^). (fm) 


3.51 


3.59 


3.67 


3.76 


^/{r'). (fm) 


3.16 


3.18 


3.20 


3.22 



> 



LU 



-6 



[440] 1/2 

[30111/2 

[312]3/2 

[303]7/2 

A 

[310]1/2 

[312]5/2 

[202]3/2 

[321]3/2 

[330] 1/2 

[200] 1/2 



Fig. 1. Neutron single-particle levels in ^"Mg labeled with the 
asymptotic quantum numbers [Nn3A\n. The solid and dotted 
lines stand for the positive and negative parities. The chemical 
potential A is indicated by the two-dotted line. 



are prolately deformed. This is consistent with the re- 
sults calculated using the Skyrme SIII interaction except 
for 34Mg [55]. The Gogny-HFB calculation suggested that 
■^^Mg is prolately deformed but soft against (3 deformation, 
and the shape coexistence in 3^'**0Mg j31j. We can see that 
the neutron skin develops as approaching the drip line; the 
difference in neutron and proton radii a/ {t''^)u ~ Vv^^ ~ 
0.35 fm in 34Mg changes to 0.54 fm in "OMg. 

In Fig. [U we show the neutron single-particle levels 
in '*'^Mg. The single-particle states are obtained by re- 
diagonalizing the self-consistent single-particle Hamilto- 
nian of Eq. ([T]). According to the present calcula- 

tion employing the SkM* Skyrme density functional and 
the mixed-type pairing density functional, *°Mg is located 
close to the neutron drip line. 



3.2 Quadrupole vibrations 

Figure [2] shows the intrinsic isoscalar quadrupole transi- 
tion strengths in neutron-rich Mg isotopes for — 0+ 
and 2+ excitations. For the K"^ = 0+ excitation, we can 
see a prominent peak possessing about 30 and 23 in Weis- 
skopf unit below the threshold energy in ^''Mg and '^''Mg. If 
we assume the strong deformation limit |42j , these intrin- 
sic isoscalar transition strengths correspond to the tran- 
sition strengths from the ground 0^ state to the 2^ state 
built on the excited K"^ — 0+ state, and those from the 
excited = 0+ state to the 2^ state built on the ground 
0|j^ state in the laboratory frame. 

On the other hand, we obtain the collective state in all 
of the isotopes for the K'^ = 2+ excitation. 

In Table [21 we summarize the ratio of the matrix ele- 
ments for neutrons and protons normalized by that of the 
neutron and proton numbers, {M^ /M.,^) / {N / Z), for the 
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Fig. 2. Intrinsic isoscalar quadrupole transition 

strengths in 34,36,38, 40]y[g j^^. ^tt ^ Q+ (yppgj. 

panel) and 2"*" (lower panel) excitations. The ar- 
rows indicate the neutron emission threshold en- 
ergies. The one-neutron emission threshold en- 
ergy is £'th,in = 3.86 MeV and 3.09 MeV, and 
the two-neutron emission threshold energy is 



-Eth,2n= 4.81 MeV and 
''"Mg, respectively. 



3.12 MeV in '^''Mg and 



Table 2. Ratios of the neutron and proton matrix elements 
Mj^/AItt divided by N/Z for the lowest excited states in Mg 
isotopes for the quadrupole if — 0^ and 2+ excitations. 







^^Mg 


3'^Mg 


3«Mg 


«Mg 




= + 

= 2+ 


1.57 
1.41 


1.58 
1.41 


1.82 
1.55 


1.91 
1.79 



lowest excitation modes. As approaching the neutron drip 
line, the contribution of the neutron excitation becomes 
large. This is one of the unique features of the excita- 
tion modes in drip-line nuclei and it is understood as fol- 
lows: In drip-line nuclei, the neutron 2qp excitations dom- 
inantly take place outside of the nuclear surface. There- 
fore, the transition strengths of the 2qp excitation of neu- 
trons become large. The proton p-h excitations, however, 
concentrate in the surface region. Consequently, coupling 
of the excitations between neutrons and protons becomes 
smaller, and the transition strengths of neutrons (A/^) and 
protons (M^) become extremely asymmetric. In Sec. 13. 2. l| 
we discuss in detail the microscopic structure of the low- 
frequency K'" = 0^ modes, and show the spatial structure 
of the excitations of neutrons and protons. 



3.2.1 Soft R-^ = 0+ modes 

In Fig. [3l we show the low-lying excitation spectra for the 
K'^ = 0+ states. Here excitation energies are evaluated 

by m 



E{I, K) = /iWRPA 



2Jtv 



(/(/+l)-i^2 



(12) 
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Fig. 3. Low-excitation energy spectra of 



34,36,38,40 



Me 



the Fermi level, one obtains a low-lying mode possessing 
enhanced strengths both for the quadrupole p-h transition 
and for the quadrupole p-p (pair) transition induced by 
the pairing fluctuations. The up-sloping and down-sloping 
orbitals have quadrupole moments with opposite signs. 

Generation mechanism of the soft — 0+ mode in 
deformed neutron-rich nuclei is understood essentially by 
the schematic two- level model in Ref. [55]. They consider 
the case where only two AA components are present in the 
wave functions both of the ground 0^ and of the excited 
states; 



0^ 



= 0+) = 



+ 62 



lAiAi 



Va2 + 62 



IA2A2) (13a) 



\K'' = 0+) = -- 



lAiAi 



IA2A 



2)- 



(13b) 



The transition matrix element for the quadrupole operator 
in terms of the vibrational frequencies wrpa and the Thoulessis then 



Valatin moment of inertia JTtv calculated microscopically 
by the QRPA as described in Ref. [44]. As we can see in 
this figure, appearance of the soft K'^ — 0+ modes is quite 
sensitive to the neutron number. 

In Ref. [44], we have discussed the generic feature of 
the low-lying K'" = 0"*" modes in deformed neutron-rich 
nuclei: In a deformed system where the up-sloping oblate- 
type and the down-sloping prolate- type orbitals exist near 



20 



0^) 



2ab 



52 



[(Ai|Q2o|Ai)-(A2|Q2o|A2)] (14) 



and it is proportional to the difference in the quadrupole 
moments of the individual orbitals composing the 0^ states. 
In the case that the quadrupole moments of the orbitals 
have opposite signs to each other, this matrix element be- 
comes large. This situation is realized in the level crossing 
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Fig. 4. Transition densities of protons (left) and neutrons 
(right) to the 7^" = 0+ states at 2.65 MeV in ^''Mg (upper) and 
at 2.30 MeV in *°Mg (lower). Solid and dotted lines indicate 
positive and negative transition densities, and the contour lines 
are plotted at intervals of 3 x 10~* fm"''. The thick solid lines 
indicate the neutron and proton half densities of the ground 
state. They are 0.054 and 0.036 fm~^ for neutrons and protons, 
respectively in ^"Mg, and 0.055 and 0.032 fm"^ in '^''Mg. 



region, where the up- and down-sloping orbitals exist. As 
the number of components increases in the QRPA calcu- 
lations, the wave function becomes more complicated. It 
is discussed in Ref. [44] . 

In what follows, we investigate the microscopic struc- 
ture of the low-lying K"^ = 0+ states, and discuss sensitiv- 
ity to the location of the Fermi level of neutrons. Because 
the deformation properties of the Mg isotopes under in- 
vestigation are not very different, Figure [T] is used for un- 
derstanding the shell structure around the Fermi level. 

In 34Mg, we obtain the collective K"^ = 0^ mode at 
2.65MeV [15]. The transition strength is enhanced by 10.6 
times as compared to the unperturbed transition strength. 
This mode is generated by many 2qp excitations, and 
among them the 2qp configurations of (j^[202]3/2)^ and 
(i^[321]3/2)^ have main contributions with weights of 0.44 
and 0.34, respectively. The chemical potential is located 
between these two levels. They are a up-sloping and a 
down-sloping orbitals, respectively. 

In '^^Mg, wc obtain two weak-collective states. The sec- 
ond = 0+ state at 3.84 MeV is mainly generated by 
the 2qp excitations of (i/[321]3/2)2 and (i^[312]5/2)2 with 
weights of 0.48 and 0.35. The lowest K"" = 0+ state at 3.17 
MeV is analogous to the collective state in ^^Mg: This is 
mainly generated by the 2qp excitations of (j/[202]3/2)^ 
and (:/[310]l/2)2 with weights of 0.58 and 0.21, which are 
a up-sloping and a down-sloping levels, respectively. 

We also obtain two weak-collective states in "^^Mg at 
3.00 MeV and 4.05 MeV. The lower state has a similar 
structure to the lowest state in -^^Mg: This is mainly gener- 
ated by the 2qp excitations of {iy[202]3/2)^ and (^[310]l/2)2 
with weights of 0.11 and 0.59. Furthermore, the 2qp ex- 
citation of (i/[312]5/2)^ has an appreciable contribution 
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Fig. 5. Transition strengths for the quadrupole-pair excitation 
in ^''Mg. The unperturbed 2qp transition strengths are shown 
in the lower panel. 



of 0.20. This excitation, however, acts destructively to the 
above excitations. Therefore the transition strength to the 
lowest state is not enhanced. The second — 0+ state 
is located just above the continuum threshold. This state 
is mainly generated by the excitations of (i^[303]7/2)^ and 
(i^[312]5/2)2 with weights of 0.38 and 0.34. 

In ■^''Mg, we can see a prominent peak at 2.30 MeV. 
This state is mainly generated by the excitations of (i^[310]l/2)^ 
and (i^[303]7/2)2 with weights of 0.50 and 0.29. The j/[303]7/2 
orbital is a up-sloping level stemming from the I/7/2 or- 
bital. Furthermore, many 2qp excitations of neutrons co- 
herently participate in generating this soft if = 0+ mode. 

The transition densities to the soft K'^ = 0+ modes 
in ^"^Mg and ''"Mg are shown in Fig.|4j Inside and around 
the nuclear surface denoted by the thick lines, neutrons 
and protons coherently oscillate along the symmetry axis 
(z— axis), indicating the P vibration. Furthermore, in neu- 
trons only, we can see a spatially extended structure. The 
enhanced transition strength of neutrons is due to this 
spatial extension of the quasiparticle wave functions. 

In order to see another interesting feature of the soft 
if = 0+ mode in ^°Mg, we show in Fig. [5] the strength 
distributions of the quadrupole-pair transition defined by 
the operator 



20 



drr^Y2o{f)i^t{r T)V'i(r i). 



(15) 



The transition strength to the lowest state is enhanced 
about 28.4 times of the transition strength of the 2qp ex- 
citation (i/[310]l/2)2. The lowest if = 0+ state in ''"Mg 
thus has a collective nature both in the p-h and p-p (pair- 
ing) channels. 

3.2.2 K"" = 2+ modes 

Next, we investigate the microscopic structure of the K"^ — 
2+ modes appearing at around 3 MeV. 
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In ^"^Mg, two quasiparticle levels near the Fermi level, 
[202] 3/2 and [321] 3/2, play a dominant role in generat- 
ing the if = 2+ mode at 3.18 MeV: 2qp excitations of 
i/[202]3/2(8)i/[220]l/2and iy[310]l/2(g)zy[321]3/2 have large 
contribution with weights of 0.51 and 0.14. 

Because a [312)5/2 level is located close to the Fermi 
level in ^^Mg, the 2qp excitation of i^[310]l/2 (g) zy[312]5/2 
becomes a dominant component with a weight of 0.54. 
This component is not large in '^^Mg, which has 5% of the 
contribution. In ^^Mg, both of the [312] 5/2 and [310] 1/2 
levels are close to the Fermi level. Therefore, the 2qp 
excitation of i/[310]l/2 (g) :/[312]5/2 has larger contribu- 
tion possessing a weight of 0.75. The 2qp excitation of 
z^[310]l/2 (g) i/[321]3/2 still have an appreciable contribu- 
tion with weights of 0.14 and 0.06 in ^^Mg and ^^Mg. 

In ^°Mg, the 2qp excitation of z^[310]l/2 (g) i/[312]5/2 
has a large contribution with a weight of 0.24 similarly to 
the K'^ = 2+ states in ^^Mg and ^^Mg. Furthermore, be- 
cause the [303] 7/2 level is close to the Fermi level, the 2qp 
excitation of ;/[312]3/2 (g z^[303]7/2 becomes a dominant 
component possessing 37% contribution. 

In all of the = 2+ states in Mg isotopes under 
investigation, the proton p-h excitations also play an im- 
portant role. Especially, the 7r[211]3/2 ^ 7r[211]l/2 ex- 
citation has an appreciable contribution with weights of 
0.21, 0.16, 0.09 and 0.09 in 34,36,38,40^g^ respectively It 
is noted that the 7-vibrational mode in ^^Mg is mainly 
generated by the neutron and proton p-h excitations of 
[211]3/2^ [211]l/2 [H]. 
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6. Intrinsic isoscalar octupole transition strengths in 
■^^■■*°Mg for = 0" (left panel), K" = 1" (middle panel) 
2~ (right panel) excitations. 



3.3 Octupole excitations 

Figure [S] shows the intrinsic isoscalar octupole transition 
strengths. For the i^^ = 3~ excitations, there are no 
peaks representing strengths greater than 1 W.u. (1 W.u. 
is about 70-90 fm^ in ■^''~^°Mg) in the low-excitation en- 
ergy region below 5 MeV. 

Because the number of 2qp negative-parity excitations 
in this mass region is small as we can see the single-particle 
energy levels in Fig. [1] none of the excited states shown 
in Fig. [H] are collective. For instance, both of the lowest 
K'^ ~ 0^ states in ^^Mg and ^^Mg are generated by the 
i^[202]3/2g)z^[321]3/2 excitation dominantly with a weight 
of 0.98. The unperturbed transition strength of this neu- 
tron 2qp excitation is about 20 fm^. About 50% of the 
transition matrix element come from the coupling to the 
giant resonances around 10 MeV and 30 MeV. 

In 3^Mg and '^''Mg, the lowest K'' = 0~ state is gener- 
ated by the i^[440]l/2g) i/[310]l/2 excitation with weights 
of 0.72 and 0.88, respectively. In ^^Mg, the iy[202]3/2 g) 
z/[321]3/2 excitation still has a small contribution of 0.08. 
The transition strengths for the K"^ = 0^ excitation in 
^®Mg and ^°Mg are about 5 times larger than in 34>36]y[g 
This is because the unperturbed transition strengths of 
the neutron 2qp excitations below 10 MeV become ex- 
tremely large due to the spatial extension of the quasipar- 
ticle wave functions around the Fermi level; the transition 
strength of the ;/[440]l/2(gi/[310]l/2 excitation is 313 fm^ 
and 720 fm^ in ^^^^^Mg. Since the i^[310]l/2 orbital is a 



particle-like level in '^^Mg, the p-h transition strength is 
smaller than in ^"Mg. About 80% of the transition ma- 
trix element come from the neutron 2qp excitations below 
10 MeV, and coupling to the giant resonance is small. In 
drip-line nuclei, the proton contribution is also small as 
in the case for the quadrupole excitations. In '^''Mg, the 
proton transition strength B{E3; 0+ ) has only 

5.0 e^fm^ whereas it has 33 e^fm^ in ^^Mg. 



4 Summary 

We have investigated the microscopic structure of the low- 
lying excitation modes systematically in neutron-rich Mag- 
nesium isotopes close to the drip line by means of the 
deformed Skyrme-QRPA method. 

In the spherical neutron-rich nuclei, the effect of the 
dynamical pairing has been investigated in detail in Refs. [45i 
H51IT7] . In Ref. f331, we showed the importance of the 
dynamical pairing in neutron-rich deformed systems. In 
a deformed system where the up- and down-sloping or- 
bitals exist near the Fermi level, one obtains the low-lying 
mode possessing extremely enhanced strengths both for 
the quadrupole p-h transition and for the quadrupole p-p 
(pair) transition induced by the pairing fluctuation. 

In this article, we demonstrated the importance of the 
dynamical pairing in deformed drip-line system. Not only 
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34 n 



Mg but in "^^Mg, the soft K'^ = 0+ modes are gener- for the K'^ = states and 



ated by the neutron 2qp excitations of the up-sloping and 
down-sloping orbitals. And we found that the transition 
strength to the soft K'^ = 0^ modes is much enhanced 
not only for the quadrupole p-h excitation but also for 
the quadrupole p-p excitation. 

We obtained the collective K"^ = 2+ excitations in 
all of the isotopes under investigation. As the neutron 
number increases, the contribution of neutron 2qp ex- 
citations to the K"^ = 2+ mode changes gradually ac- 
cording to the location of the Fermi level of neutrons. 
Furthermore, we found that the proton p-h excitation of 
[211]3/2 [211]l/2 creating the 7— vibrational mode in 
^''Mg keeps playing an important role in generating the 
collective = 2+ modes even in the near-drip-line nu- 
clei. 

For the octupole excitations, we could not find any col- 
lective modes in the low-energy region. Despite the weak 
collectivity, the transition strengths are enhanced due to 
the microscopic mechanism; In ■''*Mg and "^^Mg, coupling 
to the giant resonances at Ihwo and Shioo gives rise to 
the enhancement of the transition strengths. In ^^Mg and 
^°Mg close to the drip line, because the quasiparticle wave 
functions near the Fermi level have the spatially extended 
structure, the transition strengths of the neutron 2qp exci- 
tations become extremely large. However, the proton con- 
tribution becomes small in near-drip-line nuclei because 
the spatial overlap between the proton p-h excitations and 
neutron 2qp excitations is small. 



= (o|[x„,63i]|o) 



21 

— / pdpdzgo{p,z){4:z'^ 



3p2) (A-4) 



for the = 1 states. Here 



= 'i^f2m^n,sm, (A-5) 



1=1 



(A-6) 
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